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Abstract 

We present a calculation of the power spectrum generated in a classically symmetry- 
breaking 0(^N) scalar field through inflationary quantum fluctuations, using the large- 
limit. The effective potential of the theory in de Sitter space is obtained from a gap 
equation which is exact at large N. Quantum fluctuations restore the 0(^N) symmetry 
in de Sitter space, but for the finite values of of interest, there is symmetry breaking 
and phase ordering after inflation, described by the classical nonlinear sigma model. The 
main difference with the usual cosmic texture scenario is that the initial conditions possess 
long range correlations, with calculable exponents. 
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1. Introduction 

There are currently two broad categories of theories of cosmological perturbation: 
those which generate structure via quantum fluctuations during inflation [f], and those 
which generate structure upon the process of symmetry breaking and field ordering through 
the Kibble mechanism [2] . Inflation is a brief phase of superluminal expansion in the early 
universe, which is invoked to solve the horizon problem, the flatness problem and other 
puzzles in the standard cosmology. The existence of an event horizon gives rise to quantum 
fluctuations in the scalar held, which can be attributed to the Hawking temperature of 
de Sitter space. These fluctuations are then stretched out to cosmological scales. As the 
wavelength of a given mode expands beyond the Hubble radius H~^, its amplitude becomes 
frozen. A broad-band spectrum of fluctuations with this amplitude is generated over the 
course of the inflation. Once inflation is over, the Hubble radius grows faster than the scale 
factor a(t). The fluctuations stretched out of the the Hubble radius during inflation re- 
enter after inflation with nearly same amplitude, giving an almost scale-invariant spectrum 
of perturbations [3], [4]. Requiring that the density perturbations are consistent with the 
observed level of structure in today's universe requires some flne-tuning of parameters in 
most inflationary models. 

The symmetry breaking category [2], includes the theories of cosmic strings [5], global 
monopoles and textures [6]. A simple class of theories with broken global symmetry is 
provided by the "0(Ar)" models, where a global 0{N) symmetry group is spontaneously 
broken to 0{N — 1) by a A'"-component scalar fleld For N = 1, one has domain walls, 
N — 2, global strings and = 3, global monopoles. For A = 4, there are textures, and for 
A > 4, "nontopological textures" . In the limit of large A, there are no topological defects, 
and the theory is exactly solvable [7]. In these pictures, the universe is usually assumed to 
begin in a hot, smooth and homogeneous state, perhaps created by inflation or just special 
initial conditions. As the universe cools, the symmetry is broken and a disordered phase 
forms, containing defects of the appropriate type. As the universe expands, the tangle 
of defects is ordered on progressively larger scales. Density perturbations of constant 
amplitude at horizon crossing are generated in the process, resulting in a scale invariant 
power spectrum of perturbations. The phase ordering process involves nonlinear physics, 
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which leads to a nongaussian pattern of fluctuation [8], [9] through which they may be 
distinguished from inflation. 

At the current stage of a rapidly developing field, none of the simplest versions of the 
above theories provides a compelling match to the observations. Several variations of the 
simplest inflationary model have been introduced in an attempt to improve matters: tilting 
the power spectrum, introducing gravity waves, adding hot dark matter or a cosmological 
constant, or considering open inflation. It is also surely worthwhile to consider other op- 
tions. One of the simplest is just to combine the topological defect and inflation pictures. 
One possibility is to have standard inflation, with defects produced after reheating in the 
usual manner. Such a theory would have inflationary and defect-induced fluctuations, in 
some linear combination, uncorrelated with each other. A detailed study of reheating after 
inflation is recently been introduced by Kofman, Linde, and Starobinsky [10], where para- 
metric resonance may rapidly transfer most of the energy of an inflaton field to the energy 
of other bosons. This process may result in copious production of topological defects. Per- 
haps more interesting is the case where exponentially large defects are generated through 
phase transitions during the last stages of inflation [4]. Vishniac, Olive and Seckel [11] 
have proposed a class of models in which the inflaton is coupled to the string-producing 
field. And the strings are formed late in infiation as the infiaton rolls towards its zero- 
temperature value. The case that cosmic strings and textures are produced in inflation via 
quantum creation or sufficiently high reheat temperature has been considered by Hodges 
and Primack [12]. A second-order phase transition during inflation is studied as a natural 
mechanism to produce topological defects by Nagasawa snd Yokoyama [13]. Finally, quan- 
tum creation of topological defects through instanton tunneling during inflation has been 
investigated by Basu, Guth, and Vilenkin [14] . 

In this work, we aim at an analytic treatment of this problem using large N techniques, 
developing earlier work by R. Davis and one of us [15]. We study a renormalisable sym- 
metry beaking 0{N) scalar field theory in the large TV limit. Wc study the phenomenon 
of symmetry restoration in de Sitter space [16], and match the field modes to those in the 
subsequent radiation dominated era. This provides a set of initial conditions for the usual 
sigma model evolution, but with long range correlations. The primordial scalar field power 
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spectrum is obtained as a function of the scalar self-coupling. 

2. Theory in de Sitter Space 

We study the theory in four dimensions with the lagrangian density 
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(1) 



2.1 Large N Approximation 

The 1/N expansion of the effective potential is constructed by summing the "super- 
daisy diagrams" [17] in the limit of A?" — > oo with fixed coupling constant A. These formulas 
typically display richer structures than the corresponding leading-order expressions in or- 
dinary perturbation theory because that the leading 1/N approximation preserves much 
more of the nonlinear structure of the exact theory than does ordinary lowest-order pertur- 
bation theory. Large N behavior of a massless scalar field in fiat spacetime was first studied 
by Coleman, Jackiw and Politzer [18]. As well known, the one-loop radiative correction 
to the effective theory dynamically breaks the 0{N) symmetry [19]. However, a tachyonic 
excitation was found in the symmetry breaking state at large N, which cast some doubts 
on the validity of the large approximation. Later, Abbott, Kang and Schnitzer [20] and 
Linde [21] pointed out that the effective potential of the scalar field is a double- valued 
function, with one branch containing the symmetry breaking phase and the other including 
the symmetric solution. It was shown that the symmetric state always had lower effective 
energy. Thus the ground state of the theory is 0(iV)-symmetric, and the spontaneous sym- 
metry breaking is not possible in the large- limit. No tachyon exists in the symmetric 
vacuum, and the excited states can be studied through perturbative methods. 

This poses a problem for us - in the large N limit, the theory we are discussing does 
not have a symmetry breaking phase, even in Minkowski space. So there would be no 
classical phase ordering following inflation. However, we shall take the attitude that we 
are not in any case really interested in very large A^, (not least because the amplitude 
of density perturbations goes to zero as l/\/iV), but we merely wish to use the large N 
approximation as a tool to approximate a theory with N = 3,4 or 5. We know that at 
these modest values of N the approximation fails in Minkowski space, so there is symmetry 



breaking. It is nevertheless reasonable to hope that the approximation is valid in de Sitter 
space, where one expects symmetry restoration to occur for any N > 1 [16]. We shall 
therefore treat the problem in the quantum mechanical large approximation in de Sitter 
spaces, but the classical large N approximation (which does have symmetry breaking) in 
the subsequent radiation era. This may seem a schizophrenic point of view, but it is one 
we are forced to in order to have a chance of reproducing even the qualitative physics 
correctly. 

2.2 de Sitter Ground State 

We work in the flat slicing of de Sitter spacetime with conformal coordinate r] and 
metric 

dr^ = a^{r])[dr]^ - {dx.f] . (2) 
where the scale factor takes a{t) ~ e^^, or 

a(^) = -;^, ^<0. (3) 

We rescale a{rj) to be unity in the beginning of inflationary epoch. Thus at the end of 
inflation, or the transition to the radiation dominated era, we have ryo = —'^transition ^ 
l/H. 

In curved spacetime, the massive scalar field with mass m can be expanded as 

r (x, ry) = (27r)-3/2 j d3/c[</,,(r7)a^e*- + rMK^ e-^''"'] , (4) 

where 

[a^,a^t]^(27rfrt3(k-q). (5) 

The mode function (pi^ {rf) can be written in a separated form as [22] {rf) = a^^~"')/^Xfc iv) ■ 
The canonical quantization conditions reduce to a condition on the Wronskian of the 
solutions Xk'- Xk^vXl ~ Xl^r^Xk — ^- usual, we user the Bunch-Davies vacuum [23] for 
the scalar field: 
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where n is the spacetime dimension. The two-point correlation function can be expressed 
in dimensional regularization [24] 



8^ +zy(n))r(2 -zy(n)) 



2.3 l/A?^ Approximation in de Sitter Space 

The summation of superdaisy diagrams in flat space time has been discussed exten- 
sively [25]. The similar calculation is carried out here in de Sitter metric. We merely state 
the gap equation here, which follows from identical diagrammatics to the flat spacetime 
case 

^^-^^W^ - Ml(<t'') = ml + i„R + + ^(^^>|„„ , (8) 

where Mj^ is the effective mass squared of the theory in de Sitter space. Note that in the 
gap equation (8) (f) is the classical held, while (0^) is the two-point function given by equa- 
tion (7). From equation (7), the two-point function diverges at n = 4. Its renormalization 
is carried out in the standard procedure through dimensional regularization and minimal 
subtraction [24]: 
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A -A„+ 614^^(2 -./2), (9) 
« = {„ + ^d(d-l)(i-0M^r(2-<i/2). 

We shall only study the simplest case: zero conformal coupling, since the extension to 
nonzero ^ is straightforward. 

An alternative approach to remove the divergence in the two-point function is to apply 
an effective cutoff on the comoving momentum due to the inflationary expansion [26] . We 
shall take this approach to verify the consistency and validity of renormalization scheme 
in our model. 



After renormalisation, the gap equation (8) reads 
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where 



g = 



A ^2^9 Ml 



967r2 + A' 4 i?2 • 

The effective potential can be computed by integrating Ml 



vm = J ^d<j> + v{o) 



= j cf> . Mlicl>^) . dcl> + ViO). 
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3. FRW Solution and Mode Matching 

In the radiation dominated era, the field ordering is described classically by the non- 
linear sigma model [7] . The evolution of the N-component scalar field obeys 

To 
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In the large N approximation, Tq may be treated as a constant once scaling behaviour is 
attained for the field ordering. In momentum space, the decoupled mode functions satisfy 

To 



d'^(f>k _^ a d(j)k _^ ,^ 
d{kri)'^ kr] d{k'r]) 



{kij)' 



= 0, 



(14) 



where a = 2^^^^^^^^^ with a = 2 for radiation-dominated era and a = 4 for matter- 
dominated era. Note that equation (14) is symmetric under rj — > —rj. In the radiation 
dominated era, the solution of (14) can be generally expressed as 
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After the end of inflation, the mode functions (6) of in de Sitter space is matched 
to the FRW classical solution (15) by requiring that scale factor 0(77), (l)kiv)i 4>kiv) 
be continuous at the transition point r] = —r]o . A simple model of the transition can be 
constructed by considering a metric of the form 



a{r]) 



{ri + 2rio)/Hril 



and the matching condition reads 
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V < ~VOj de Sitter 
V> -Vo, FRW 
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(16) 



(17) 



where dot denotes d/dijQ. The solution of (17) is expressed as 



Ci = 

C2 = 

where the coefficients read 
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v'oHi^Hkiio)B - Hl^\kvo)C 
rjlHl^\krio)A - Hl^\krio)C 



(18) 



^ = iy(i)(/cr7o)-— iy«(/cr/o), 
Zrjo 

B = Hi^\kno) - :^H(^\kvo) , 
C^v'oHi'\kvo) + lvoHi'\kvo) 



(19) 



Easily we can verify the normalization condition for Bunch-Davis vacuum. The normal- 
ization conditions for the FRW mode functions are also automatically satisfied. 

The mode function (15) in FRW can be separated into a 'classical' part (by which 
we mean the component which survives with large amplitude long after the transition to 
radiation domination) and a 'quantum' part (which diminishes quickly after infiation): 

(f>k{'n) = (f>cl + (Pq , 

= a{k) J^{kr]) + b{k) Y^{kr]) (20) 



where the coefficients a and b write 
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At the epoch of most importance for structure formation, or when a given mode re-enters 
the Hubble radius, krj ^ 1, the 'classical' part dominates the 'quantum' part by a factor 

2/i 

with fi > 0. Consider how small rjo is, this is an extremely large number. It 
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is therefore reasonable to neglect the quantum part of the wavefunction at FRW space. 
Therefore, the classical non-linear sigma model is a valid approximation to the 0{N) scalar 
field at FRW epoch. 

After matching the quantum fluctuations of the scalar field to the classical evolution 
of the nonlinear sigma model, we can now relate the original parameters in the quantum 
field theory to the value of 0^ = 0q in the classical theory. We have 



/J, = 2 — v . 



(22) 



And at large rj ^ rjo 



d^k 1 



(27r)3?7^^' ' 

r2(z.)r2(^)r(/i-i/2) 



(23) 



47r4 T{ii + 3/2) 

At later time in FRW, the classical dynamics of the theory is totally fixed by the tree 
diagrams. Thus, the renormalised values and A are fixed by 

A 



m 



6N 



(24) 



where 4>o{H) is a function of H as fixed by (23). To re-emphasize the point, we are treating 
the theory as classical in the post-infiation era, with the given renormalised parameters, 
with the justification that for modest values of N and A this should be a reasonable 
approximation. 



4. Eflfective Potential and Vacuum in de Sitter Spacetime 
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The effective potential of the theory in de Sitter space is studied under the as- 
sumption of equation (24). After plugging (24) into the gap equation (10), our theory- 
will have only one free parameter — the renormalised coupling constant A. The effective 
potential Ve// is obtained for different regions of A and the corresponding vacuum states 
are found by minimizing Ve// 

d(f> ^ del)'' (25) 

= cl>Ml{ct>'). 

To satisfy equation (25), the vacuum state will either have |0| =0, where the symmetry is 
restored or Mj^ = 0, where the symmetry is broken. 

4.1 Weak coupling: X <^ N 

Weak coupling implies the small mass limit: \M^J ^ H^. At first, let's study the 
weak interacting limit of the gap equation (10) with X ^ N and jM^^I <^ H^. The leading 
order terms in this limit 



, _ ,/9 Ml _ 3 IMl Ml 



00 = 



3NH^ H' 
lo^Ml 



o A ,0 XH H 



(26) 



^N^^ 327r2 Ml 



2 ^'ZH-^' 2 Ml 

The leading order of the gap equation (10) is 



Ml _ I X 



i?2 V 327r2 ■ 
The variable v{X) can be expressed as 

3 1 IX 



(27) 



2 127r V 2 

Generally, the relation of Ml{(jp') over (jp' is plotted at Fig.(l). At the region Ml < A, 
(jp' is negative, which is obviously unphysical. Thus, Ml cannot be zero and the ground 
state of the scalar field is the 0(A'")-symmetric state 

101 =0, Ml = A>^. (28) 
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The effective potential is obtained by integrating M^, as sliown qualitatively in Fig. (2). 
For small 0, Ve/ / can be easily calculated 

VefficP) = I 4>Mld4>^ j ci>{A + l^</>2)d</,, 

(29) 

= 2*^ + m^^XN^^ ■ 

For (j) > (j)max, the effective mass squared Mj^ turns complex and the effective potential 
has an imaginary part, which is related to the particle creation in the scalar field cj) [21]. 
The imaginary Vef / and the vacuum instability imply that one can actually deal with the 
vacuum state with the field only during such a short time interval, for which the energy 
density of the particles created during this time, will be much less than Re(Ve//), and the 
effects of the tachyonic instability of vacuum can be neglected. 

The primordial power spectrum of the scalar field is obtained by taking the long- 
wavelength limit of the two point correlation function, 

P(/c,r/)= mriW-M) , 

= a{kfJl{kr^), 

1 (30) 

4.2 Intermediate Coupling: X N 

As A is increased to the order of of N, the effective mass squared turns to be a 
double- valued function of (jp, as illustrated in Fig. (3). A new 0(Ar)-symmetric solution C 
emerges with Mj^ ~ H'^. Vacuum solution C has a lower effective potential compared to 
vacuum A. This can be shown by taking the integration of from A to the transition 

point B and then from B to C 

B c 
ViC) = J <t>Ml{4>)d4> + j 4>Ml{4>)d4> + V{A) < V{A) . (31) 

A B 

Thus, state C is the true ground state, as illustrated in Fig. (4). The large effective mass 
vacuum state is an unique solution in the large N approximation, which is absent in the 
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traditional one-loop analysis [27]. At the limit of — > |, we obtain almost Gaussian 
white noise spectrum with z/ — > 



Ml 9 2 2 1 



II ' 4 ' 907r4zy2 " 

The leading order of the gap equation (10) is 

2. , , , 4(967r2 + A) ^ .oo^ 
zv"' logiu) + — ^ = . (32) 

^ 87-907r4 

Numerical simulations show that small ly solution is obtained around X N and the 
vacuum state C exists at the range oi < v < 0.8. And the power spectrum is 

1 1 



J^2i' J^O ~ 1.6 

For very large A, or A ^ N, the contribution from high order terms will be significant 
and the the 1/N expansion is no longer reliable [21]. The power spectrum index 2^ versus 
A is plotted in Fig. (5). 

5. Verification of the Renormalization and Matching Schemes 

The divergence of two-point function (0^) in de Sitter space can be removed by a less 
formal but more intuitive methods: inflationary stage naturally provides both ultraviolet 
and infrared cutoffs at momentum space. Let's consider the physical wavelength of mode 
k, Xphy = a{r])/k. Any mode with physical wavelength greater than the Hubble radius at 
the beginning of inflation {rj = rji), or 

Ap/ij/s > 1/H ^ k < H, with?7j = and a{r]i) = 1 , (33) 

11 

would be totally stretched out by inflation and is irrelevant to the physical processes after 
inflation. Also, any mode function with physical wavelength shorter than the Hubble 
radius after inflation will get more and more subhorizon in FRW, so it does not contribute 
to (0^) either, 

Khys < i/H ^ k > 1/r/o . (34) 

The same scheme was applied by Vilenkin [26] to study the quantum fluctuation of a 
scalar fleld during inflation. 
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We are most interested in the small coupling region, where <^ 1 and ~ § — | 
From equation (6) 

l/»?o 1/t?o o 



(27r)3 

NH^V\v)4?^ } (35) 
~ ;;^2-^ / d{kr]o){kvoy 



7r2 3-2z/ ' 

In FRW space, the mode functions satisfy equation (23) 

^2 ,, 2ifV ,^,^_,^, _,,, r(//-i/2) 
<^o-A^^^(r(-)r(/.)) r(^ + 3/2) ' 

9Arii'2r2(zy) 1 



(36) 



The (0^) terms match remarkably well from de Sitter stage to FRW stage in this limit. 

Around the Gaussian white noise limit, we have ~ 9/4, ~ 0, and (l)v{x) ~ 
^ln{^). The de Sitter solution has 

i/^ 1 

d^k ,o,, . NH^ r .o, 2 



^^'^ = / / ^(^^o)(/c.o)^^n^(^) . (37) 

H 

And the mode functions of FRW satisfy equation (23) 

^S-^-^/i^'-^li^)^.^!**). ^^^^ 

y d{kvo){kvofln\-^) . 
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Equation (37) and equation (38) are dominated by small argument integrand, so they 
match each other relatively well. 

The consistency among above two schemes of renormalisation verifies the soundness 
of our approaches. 

6. Conclusion 
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In conclusion of the paper, topological defect originated from quantum fluctuations 
during inflation is a viable and promising picture for cosmic structure formation. We 
show that the de Sitter vacuum state provides a natural initial condition for the symmetry 
breaking scalar field in FRW, which is described by the classical nonlinear sigma model. 
Due to an early stage of inflation, the scalar field possesses superhorizon correlations. 
Thus, this scenario might have more large scale power than usual cosmic defects models. 
By proper adjustment of the only free parameter: the coupling constant A, we can obtain 
power spectrum ( = k~'^^) ranging from exact scale-invariant to Gaussian white noise. 
In the weak coupling case (A <^ A^), we have 2z/ = 3—5, where 5 ~ For intermediate 

coupling region {X N), 2v takes the value of to 1.6. 
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Figure Captions 

Figure 1: The effective mass squared M^^ versus in tlie iimit of weak coupling witli X ^ N. 
At the symmetric phase point A: 0^ = and M^g{A) > 0. 

Figure 2: The effective potential Veff in the weak-couphng hmit. The effective potential is 
calculated by integrating M^^. The O(A^) -symmetric ground state lies in 0^ = 0. 
For 0^ > (fyfnaxj effective potential turns imaginary representing the generation of 
classical field. 

Figure 3: The effective mass Mj^ versus (j) in the intermediate interaction range with A ~ A^. 

New vacuum solution C emerges in the limit of large effective mass. The positive |0p 
region is divided into two separated regions, from A to B, and from B to C. 

Figure 4: The effective potential Vg// in the intermediate coupling region. The 0(A^)-symmetric 
ground state is the large mass solution C which has lower effective potential. 

Figure 5: Power spectrum index 2u with respect to A. {Pt = k^'^'^) At A ^ A^, we have the 
almost scale-invariant spectrum 2^ ~ 3. The vacuum state jumps from A to C when 
A reachs the order of N. Gaussian white noise limit is obtained around X N. There 
is no solution ior X^ N. 
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